Filtering property of myelinated internode can change neural information representability and might trigger a compensatory action during demyelination

In this paper, for the first time, we showed that an Internode Segment (INS) of a myelinated axon acts as a lowpass filter, and its filter characteristics depend on the number of myelin turns. Consequently, we showed how the representability of a neural signal could be altered with myelin loss in pathological conditions involving demyelinating diseases. Contrary to the traditionally held viewpoint that myelin geometry of an INS is optimised for maximising Conduction Velocity (CV) of Action Potential (AP), our theory provides an alternative viewpoint that myelin geometry of an INS is optimised for maximizing representability of the stimuli a fibre is meant to carry. Subsequently, we show that this new viewpoint could explain hitherto unexplained experimentally observed phenomena such as, shortening of INS length during demyelination and remyelination, and non-uniform distribution of INS in the central nervous system fibres and associated changes in diameter of nodes of ranvier along an axon. Finally, our theory indicates that a compensatory action could take place during demyelination up to a certain number of loss of myelin turns to preserve the neural signal representability by simultaneous linear scaling of the length of an INS and the inner radius of the fibre.


Results
INS behaves as a low-pass filter with f L higher than biologically relevant range in normal oper- ating condition Equation (12) (see Section "Materials and methods") implies that an INS is a 2nd order system, and its frequency response should be of bandpass type, i.e., it will pass a certain frequency range unattenuated while attenuating the other frequencies that do not fall within that band.A straightforward way to determine such a passband is to examine the pole-zero plot of the system Transfer Function (TF; Section "Materials and methods") for each of the fibres as shown in Fig. 1a.As expected, each TF results in two poles and one zero in accordance with a bandpass system's characteristic.But, surprisingly, for all the fibres, their only zeros coincide with one of their two poles compensating each other's effect leaving only the other pole to dominate the system's behaviour.This effectively makes an INS a lowpass filter.More interestingly, the dominant poles of all the PNS fibres are located very close to each other implying that despite having different values of L , r and M (Table 2, Section "Materials and methods"), all the PNS fibres would exhibit similar time and frequency response characteristics.On the other hand, the dominant poles of the CNS fibres are separated significantly to show that each of them may exhibit their own characteristic time and frequency response property.To verify this and to find the f L (following the process outlined in Section "Materials and methods") for each fibre we derived the Bode Magnitude plot from each fibre's TF as shown in Fig. 1b.Consistent with Fig. 1a, all the fibres show lowpass filter characteristics with all the PNS fibres exhibiting similar frequency-dependent gain characteristics with their f L coincided.Overall, this is surprising since the individual PNS fibres have different INS geometries (see Table 2), and thus are expected to show different gain characteristics.The only possible scenario under which such behaviour can occur is if L , r and M are interrelated so that multiple but appropriate combinations of them can lead to the same f L .It is worth noting that the numerical values of g and γ (see Table 2) for all the PNS fibres are the same despite having different geometries and therefore strongly suggests such interdependency.
On the other hand, the CNS fibres show their individual characteristics with lower gain and f L compared to the PNS fibres.Interestingly, although their g-values are close to those of the PNS fibres, their γ-values differ significantly not only compared to the PNS fibres but also among themselves, which could be suggestive to their location-specific unique characteristics.However, such a different behaviour precludes a definitive conclusion that L , r and M are interrelated.But, in all the cases, f L > 1 KHz implies that no AP disruption could take place within an INS in their normal operating condition.

Representability of temporal variation of a stimulus could change depending upon the myelination state of a fibre
Next, we progressively changed M for each fibre while keeping their L and r constant and extracted f L , (Section "Materials and methods).The results are shown in Fig. 2 for the case of A α11 fibre (see Section "Materials and methods" for its definition) as an illustrative example (for the corresponding figures for all the other fibres see Fig. S1 in SI).As M decreases, f L moves toward the lower frequency band-reducing the fibre's NCC-eventually moving below 1 KHz, in this case, when M < 50.
The implication could be explained by taking the specific example of f L = 767.8Hz when M = 30 .At this condition only those states of a stimulus that could be represented by AP firing rate ≤ 767.8 Hz will pass through the INS while any state representation requiring AP firing frequency above that will be blocked within the INS-effectively filtering out part of the original neural code.An illustrative example has been shown in Fig. S2 (SI) with AP train of 1 KHz frequency.We observed the same pattern for all the fibres: f L changes linearly with M (see Fig. S3, SI) implying that to reach a certain NCC, a fibre needs a certain minimum M .But the rate of change of f L with respect to M appears to be steeper for the fibres with smaller r and L .This suggests that the locus of f L might not only be dependent on M but possibly on r and L too.Therefore, it is more prudent to study the behavior of f L as a combined effect of these three parameters.

Figure 1.
The frequency characteristics of the INS for the fibres under consideration (see Section "Geometrical parameters" in Section "Materials and methods" for explanation of the legends and Table 2 for their dimensions).(a) shows the pole-zero plots with the circles representing zeros and " × " representing the poles.Note the coincident poles and zeros for all the fibres close to the origin.(b) The frequency dependent gain characteristics of all the fibres.Filled circles denote the f L for each fibre calculated following Eq.( 13), see Section "Materials and methods" Section "INS system modelling and f L calculation".All the PNS fibres exhibit identical characteristics and their f L is coincident while the CNS fibres show their individual characteristics.
At the same time, we found that the reduction of M will not affect the shape and amplitude of an AP since most of its spectral power is concentrated within a very narrow low-frequency band of 1-5 Hz with the maximum at 1.9 Hz (see Fig. S4, SI).Therefore, the shape and amplitude of an AP could only be changed significantly for f L ≤ 1.9 Hz which is at least an order lower than the lowest value of f L corresponding to M = 1.

L, r and M are linearly interdependent irrespective of the type of nerve fibres
To study how f L varies as a combined effect of r , L and M , we plotted f L for each M for all the fibres against their corresponding g and γ as shown in Fig. 3. Surprisingly, instead of being scattered around on a g-γ surface, the loci of f L for all the PNS fibres overlap with each other forming a single curve as shown by the brown dashed trajectory in Fig. 3, that could be represented by the fitted linear Eq. (1) (95% confidence bounds, R 2 = 0.99):   This means, irrespective of fibre types, f L is a linear function of g and γ , and every possible f L of all the PNS fibres could be specified by a unique tuple (g, γ ) whilst in the case of CNS fibres the same is true but in fibre- specific manner.Therefore, for a specific f L the following two conditions hold: where C 1 and C 2 are two unique constants specific to a f L for a fiber.Substituting Eqs.(7k) and (7l) (Section "Materials and methods") in Eqs.(2a) and (2b) we get Substituting r from Eq. (3a) into Eq.(3b) and solving for L we get Since all the parameters on the right-hand side of both equations are positive, the second solution does not exist since L cannot be negative, leaving only the first solution.For that solution to exist the following condition must be satisfied Since both g and γ are always < 1 and C 2 ∝ 1/M 2 , Eq. ( 5) is always satisfied for M > 0 which means that L and M are linearly proportional.Combining this with Eq. (3a), proves that r, L and M are linearly dependent on each other, and this is valid both for a PNS and CNS fibre.That means, depending on a 'desired' NCC, irrespective of the type of a fibre, r, L and M are chosen to define the geometric structure of an INS following a linear interdependency rule, that allows unperturbed AP generation at the next NR for all operating frequencies up to f L , which defines the upper limit of the NCC.
A consequence of Eq. ( 4) is, to ensure the full range representability of temporal variation of a stimulus, i.e., f L ≥ 1 KHz, there exists a limiting L/M value for each fibre which could be calculated as follows: Let (g 1 , γ 1 ) be the tuple at which f L = 1 kHz for a fibre and the corresponding constants in Eqs.(3a) and (3b) are (C 11 , C 12 ) .Therefore, from Eq. ( 4) we get, From Fig. 3, since all the PNS fibres intersect the f L = 1 KHz plane at a unique point (g 1 = 0.961, γ 1 = 0.0052) , from Eq. ( 6) L/M ≤ 0.5µ m is the condition for a PNS fibre to keep its f L ≥ 1 KHz.Similar limiting conditions for the CNS fibres could also be derived but since they exhibit their individual characteristics, such limiting condition needs to be calculated in fibre-specific way.

Loss of myelin turns under pathological conditions might trigger a compensatory mechanism to preserve a fibre's NCC by dynamically scaling r and L linearly
Since both r and L are linearly dependent on M , with reduction of M during pathological conditions related to demyelination, there must be different possible combinations of r and L for which a unique tuple (g, γ ) and hence a unique f L could remain constant if a fibre 'intends' to maintain that f L .To illustrate the scenario, let's assume that the Aα11 fibre intends to preserve its normal operating f L = 10 KHz (approx.)during such a scenario.From the values of (g, γ ) in Table 2 corresponding to f L = 10 KHz, using Eqs.(2a) and (2b), we calculate C 1 = 0.005 and C 2 = 102.041 .Applying these values in Eqs. ( 4) and (3a), we calculated the values of L and r respectively for each value of M from 400 to 1 (Fig. S5, SI).Subsequently, we derived the TFs of the INS for each of these (r, L, M) combinations and plotted the corresponding Bode Magnitude plot in Fig. 4 which shows that even in this case the f L remains constant.We have observed the same behavior for all the fibres with different f L values.The experimental findings that during demyelination the INS length is shortened when interpreted based on the results from Fig. 4 clearly suggest an on-going compensatory mechanism to maintain a certain constant f L , by simultaneous linear scaling of r and L with M following Eqs.(3a) and ( 4).Although, as per our knowledge, Vol:.( 1234567890) www.nature.com/scientificreports/no such direct experimental result is available with respect to r , our result suggests that the experimentally measured overall thinning of the outer diameter of a fibre in such a scenario must have the effect of reduction of r embedded in it.

CD remains constant during a compensatory scenario, but CV is reduced
Next, we calculate CV using the procedure described in Section "Materials and methods" by first calculating the CD. Figure 5 shows the result for Aα11 fibre as an illustration with its corresponding geometrical parameters given in Table 2. Clearly, CD falls exponentially with increasing AP train frequency.Consequently, given a fixed L , the CV is also frequency dependent and increases exponentially with AP firing frequency.But the situation is different during the compensatory scenario described in Section "Loss of myelin turns under pathological conditions might trigger a compensatory mechanism to preserve a fibre's NCC by dynamically scaling r and L linearly".Taking the same example of Aα11 fibre, but this time with M = 50 and its correspond- ing L = 0.25 mm (Fig. S5, SI), we illustrate the phenomenon in Fig. 6.Here it is evident that at the intended f L = 10 KHz, the CD remains constant at 5.371µs for both M = 50 and 400 but the CV changed from 372m/s to 44.68m/s since in this case L is reduced from 2 mm to 0.25 mm to keep CD constant.Therefore, it is not surprising that reduction of CV has been observed experimentally during loss of myelin but its interpretation that the AP signal slows down within the INS is not necessarily true-in fact it is suggestive to the ongoing compensatory action to preserve the NCC.

Discussion
Our findings that there exist multiple combinations of r, L and M for preserving the unique tuple (g, γ ) and hence, a unique f L on one hand, support the general result derived in 1 and is instrumental in qualitatively explaining several experimental observations that were in apparent contradiction with conclusions in 1 as follows.
In CNS the INS length is progressively reduced while the diameter of the NR increases progressively along the axon, 10,11,19 Since our theory suggests that multiple combinations of L, r and M could satisfy the same f L , it is not surprising to have different INS geometry along the length of an axon for preserving a certain f L .As suggested in 19 , such distribution is possibly guided by the availability of local resources.To explain the progressive increase of NR diameter with shortening INS let's assume that the typical length and radius of an INS and NR are ( L INS , r INS ) and ( L NR1 , r NR1 ) respectively.As L INS reduces, L NR1 becomes proportionally longer to L NR2 and r NR1 changes to r NR2 (say).Such a change in geometry will increase the delay in the NR from t NR1 to t NR2 while according to the findings in Section "CD remains constant during a compensatory scenario, but CV is reduced", the delay introduced at INS ( t INS ) remains constant.Therefore, considering the INS-NR as a combined unit, the total delay changes from t D1 = t INS + t NR1 to t D2 = t INS + t NR2 .Since NR is an unmyelinated section, the delay per unit length is similar to that of an unmyelinated fibre and thus, following the results of 1,20 , t NR1 ∝ L NR1 / √ r NR1 and t NR2 ∝ L NR2 / √ r NR2 .If the final goal is to keep the NCC of the fibre constant, then the condition t D1 = t D2 must hold true which implies r NR2 = r NR1 (L NR2 /L NR1 ) 2 .Since (L NR2 /L NR1 ) > 1 , r NR2 > r NR1 which qualitatively matches the experimental findings in 11 .On the other hand, if the target is to keep CV constant instead of the delay in NR, CV in an NR being proportional to √ r NR , we get r NR2 = r NR1 (t D1 /t D2 ) 2 .Since in this case (t D1 /t D2 ) < 1 , we get r NR2 < r NR1 which is the opposite to the experimental observation in 11 .

Shortening of INS during demyelination and remyelination
While our results directly explain why one might expect shortening of INS in the scenario of myelin loss ("Loss of myelin turns under pathological conditions might trigger a compensatory mechanism to preserve a fibre's NCC by dynamically scaling r and L linearly"), this also gives a possible explanation to why the INS remains short during remyelination.In 13,14 , the authors found that the INS geometry changes with learning and might be a life-long process.This suggests that during remyelination one must not expect the regenerated myelin structure like a normal population immediately as the biological adjustments occur progressively over time influenced by learning-related plasticity property.Since our results show that multiple combinations of L, r and M are possible to attain a specific f L , it is not hard to postulate that during remyelination the myelin geometry is configured only to support the 'baseline' NCC and over a longer observation window, one is likely to find more changes as suggested in 13,14 .

CNS fibres might have unique properties depending upon their location
Experimental observations in 19,21 showed that there exists regional difference in Oligodendrocyte properties in CNS.In our results described in Section "L, r and M are linearly interdependent irrespective of the type of nerve fibres", this is reflected by the similar functional form of f L for both the CNS fibres while differing by only the fitting constants which in essence encapsulate the values of the biological constants such as membrane resistivity, specific capacitance, and membrane thickness.Therefore, our theory is consistent with the above-mentioned experimental observations.On the other hand, the almost similar characteristics shown by the PNS fibres in this regard is possibly since each PNS fibre is myelinated by a single type of Schwann cell with similar property.www.nature.com/scientificreports/

Limit of compensatory action
The results presented in Section "Results" implies that in myelin loss-induced compensatory scenario it is possible to maintain the NCC by scaling down the INS geometrical parameters even when M reaches to 1 .This is counterintuitive given the fact that partial myelin loss also can lead to pathological signal propagation conditions.Therefore, there must be a limit of compensation which the theory didn't spell out explicitly.Numerical evaluation of Eqs.(3a) and ( 4) shows that given a certain f L , to sustain such compensation, r becomes < 1µm after a certain reduction of M (Fig. S6, SI).It is a well-known fact that for the PNS fibres with such a small radius the advantage of myelination is lost due to energetics, speed, and mechanical stability requirements 1,20 .Therefore, this is an obvious point to stop any further compensation.On the other hand, any compensatory action will also elongate the adjacent NR and to mitigate the associated increase in delay, it's diameter must be increased.This will bring the available space and resource requirements into question as described in 20 .These factors can certainly play a role in restricting the compensation action even earlier.However, appropriate directed experiments are needed to explore this area.

Implications on neural information processing
The traditional view of neural information processing is based on the concept of spatial neural processing-the time synchrony between AP spikes arriving through different paths in a neural circuit 13 .This time-synchrony is influenced by the INS lengths along the circuit paths as they change the CV.Our results showed that despite the observed reduction of CV, the actual CD can remain constant within the limit of compensation.Since CD is the actual parameter determining the time synchrony between different paths in a neural circuit, the abovementioned traditional view is valid only after the limit of compensation has been reached.More importantly, beyond the compensation limit, not only the time synchrony but also the representation of the actual neural information is changed due to low-pass filtering effect of an INS.This is a component that was previously ignored in developing the theory of neural information processing, particularly in the pathological conditions related ( ω L = 62832rad/s ) remains constant, in both cases of M = 400 and 50 , the CV (Bottom) in the latter is much lower compared to the former.If measured within the range of 1 KHz ( ω ≤ 6283rad/s ), the CV will be even lower and consistent with experimentally measured range although this leads to a misinterpretation of the phenomenon.
to myelin loss.Thus, our theory shows that this component must be added in developing such a theory while adjustments in spatial neural processing concept should be made in light of CD remaining constant within the compensation limit.

Effect of radial components
Our analyses were based on the DC model 18 that didn't explicitly include the effects of myelin radial components.Gow and Devaux 22 showed that its absence decreases CV particularly in small fibres that are the largest constituents of CNS.Although initially thought that these radial components have adhesive properties to hold the myelin lamellae together, electrophysiological modelling in 23 showed that they have insulative properties by formation of diffusion barriers and thereby regulating current flow through the myelin sheath.Using rigorous biophysical and ultrastructural techniques including X-ray diffraction, 24 conclusively established this diffusion barrier property of the radial components.The main reason for not including it explicitly in our analyses is the lack of appropriate data.The effect of radial components on CV for small fibres described in 22 used computational modelling techniques with two alternative models, one with the resistance resulting from the Tight Junction (TJ) of the radial component in series with myelin resistance-capacitance assembly; termed as the TJM, and the other with the resistance of TJ in parallel with myelin resistance in the conventional DC model-for ease of discussion we term this model as DCTJM.Electrical property of epithelial cell was used to characterise the electrical property of TJ.However, the works described in 19 , 21 show that in CNS the electrical property of INS is strongly dependent on local Oligodendrocyte properties.Therefore, for CNS fibres, it seems more pragmatic to characterise the electrical properties of TJ in fibre-specific manner.Currently, to our knowledge, such data is unavailable for CC and CB fibres-the two small fibres explored in this work.Despite the lack of data, following the approach described in 22 , we have explored the possible implications of radial components on variations of f L in small fibres.The native configuration of the fibres is defined by their respective geometrical parameters given in Table 2, and thickness and resistivity of intramyelinic space as d = t m and ρ = 1.96�m 23 respectively where t m is the thickness of a lipid layer (see Fig. 7, Section "Materials and methods).The main results could be summarised as follows: 1. From a system point of view, the TJM model proposed in 22 introduced an extra zero in the TF of all fibres making a change in shape of their Bode magnitude plots as shown by the green curves in Fig. S7 (SI).For illustration purpose we have only depicted the response of the largest fibre Aα11 and the two smallest fibres CC and CB.However, the location of that extra zero is far away from their corresponding f L and therefore www.nature.com/scientificreports/practically has no effect in the system's response up to their respective f L .On the other hand, the number of zeros and poles of the DCTJM is the same as the DC model used here.Effectively, all the three models imply that an INS is a lowpass filter.2. Consistent with the results of 22 we found that for large fibres ( Aα11 here) there is no effect of TJ on f L (Fig. S8, SI) in its native configuration for all the three models.For small fibres (CB and CC) the TJM and the DC model used here gives no difference in f L (Fig. S9, SI).However, in the DCTJM the f L is reduced by 22.32 Hz for CB and by 24.67 Hz for CC compared to the other two models.3. The above-mentioned phenomenon for CC and CB is strongly dependent on the value of d .For d < t m , no difference between the three models could be found (Fig. S10, SI).On the other hand, for d ≥ t m , although the TJM and the DC model show the same characteristics, marked difference has been observed in the case of the DCTJM (Fig. S11, SI).In this case f L is reduced by 150.16 Hz for CB and by 186 Hz for CC. 4. The linear interdependency of L , r and M , irrespective of the type of nerve fibres still holds true with the TJM and DCTJM model (Fig. S12, SI).
Effectively we found that even for small fibres, the TJM and the DC model used here always give the same results as far as f L is concerned but the DCTJM gives different results depending on the chosen thickness of intramyelinic space.Consistent with 18 , these results underscore the important role played by the periaxonal resistance R p .In all the three models, R p is in parallel with their respective myelin-TJ assembly resistance and therefore the systems' response is strongly dependent on their relative values.This is also reflected by the dependence of the system's response on d which determines the value of the resistance of TJ.Given the results in 19 , 21 showing the region-specific difference in Oligodendrocyte properties, and hence the fibres' electrical characteristics in CNS, these results show that fibre-specific characterization of R p and myelin-TJ assembly resistance (in essence, the value of d ) is necessary to fully characterize the effects of radial components on frequency domain proper- ties of INS in these fibres.

Limitations of the proposed theory
The first limitation of our theory is that here we considered INS as an isolated system and "clean" turn-by-turn myelin loss.In practice, loss of myelin is "patchy".This limitation could be overcome by partitioning an INS into a "normal" and a "patchy" region and then deriving TF for each region followed by multiplication of these TFs to create a composite TF.Subsequently, this TF could be used to obtain the results corresponding to those described in Section "Results".
The second limitation is, although more than 1 KHz AP firing rate has not been observed in practice, still our theory predicts that the NCC of a healthy fibre could be as high as 10 KHz.Given the observation in Fig. S6 (SI) which shows that according to the condition L/M ≤ 0.5µ m presented in Section "L, r and M are linearly interde- pendent irrespective of the type of nerve fibres", no PNS fibre could maintain f L ≥ 1 KHz for M ≤ 3 , apparently this is to increase the resilience of a fiber against possible myelin loss.However, we would like to point out that this value is dependent on the choice of R a , R m , R my , C m and C my (see Section "Materials and methods") which in www.nature.com/scientificreports/turn depend on the biological constants ρ a , ρ m , ǫ m and t m .As discussed in the Section "Materials and methods", there are considerable variations in numerical values of these parameters reported over the years.Using their different values results in different numerical values of f L .We are unable to comment on such variations of these physical constants measured in different experiments, but this wide range of variation certainly indicates that not only the f L but all the other quantitative results derived in Section "Results", can only be considered accurate with respect to the numerical values we have used in Table 1-Section "Materials and methods".However, even using different values of the above-mentioned constants the phenomena we described here remain the same.Therefore, we tried to avoid providing quantitative results wherever possible and preferred describing the essential parts of our theory in its functional form.We suggest that more experiments are needed to measure these values accurately to confirm the numerical values of the results we presented in Section "Results".The third limitation is, the compensatory phenomenon predicted by our theory, although widely speculated by several authors, is yet to be detected experimentally.Therefore, we suggest directed experiments based on our theoretical arguments to either prove or disprove such phenomenon which will give more insights into the operation of our nervous system under pathological conditions.

The modelling circuit structure and electrical parameters estimation
The structure of the INS we used for modelling purpose is shown in Fig. 7 along with the notations we used to specify its geometry.
We approach the modelling considering cylindrical symmetry at the granularity of number of myelin turns rather than the traditional approach of using outer diameter since we are more interested to understand how the system behaves during changes in myelin turns.
We used the Double Cable (DC) model 18 -which includes the effect of periaxonal space-as the main circuit modelling structure of an INS (Fig. 7).To estimate the electrical components, we assume that all the myelin lamellae are already compacted with each lamella composed of two tightly apposed lipid bilayer (axolemma membrane)-a standard approach for modelling an INS.Consequently, if each lipid bilayer has a thickness of t m , then the thickness of each myelin lamella is 2t m , and its capacitance and resistances could be expressed in terms of the corresponding parameters of axolemma as follows: where R m and C m represent the resistance and capacitance of the axolemma which, for INS length L could be expressed as where ρ m and ǫ m are the resistivity and permittivity of the axolemma respectively.On the other hand, the axo- plasmic resistance could be given by where ρ a is the resistivity of the axoplasm.Following the cylindrical symmetry, assuming that the thickness of the periaxonal and paranodal space being d 1 and d 2 respectively with their lengths being L 1 and L 2 respectively, the corresponding resistances ( R pa and R pn respectively) could be calculated as where ρ 1 and ρ 2 are the resistivity of the periaxonal and paranodal fluid respectively.The total periaxonal resist- ance is the sum of the above two resistances and thus could be expressed as The choice of the fundamental electrical constants used in this work are shown in Table 1.We observed that the experimentally reported values for these parameters are different in different published papers.A specific example is ρ a which has been reported in the range of 3.3-0.7 Ωm [25][26][27] making it difficult to determine which one to use.Therefore, in such cases we have used the mean value of the range.
considering the membrane resting potential and the membrane potential at which Na channel activates are −70 mV and −55 mV respectively, with the peak of an AP being +40 mV, to trigger an AP at an NR its membrane potential must be at least +15 mV at the peak of an AP.This leads to the following condition for triggering an AP at the NR: Accordingly, ω L = 2πf L represents the frequency at which G(ω L ) = −8.519dB ,where ω L and f L are the angular and linear frequency respectively.This could be detected by examining the Bode magnitude plots which describes the characteristics of G(ω) vs ω of a TF.

The simulation flow
Our main simulation experimental process could be described as follows-first we calculate the TF for each of the fibers using their geometrical and electrical parameters given in Tables 2 and 1 respectively and then replacing them in the circuit structure shown in Fig. 7. Subsequently, we calculated G(ω) for each of the fibres using Bode Magnitude plot and extracted the corresponding f L following Eq.( 13).This step gives the system behaviour of an INS of that fibre in its normal operating conditions.Next, we progressively reduced M for each fibre from their values given in Table 2 down to 1, which represents the limit of myelin loss to simulate the scenario of changes in myelin turns.At each step of this process, we repeat the process of calculating TF and corresponding f L that shows how f L changes with respect to the changes in myelin turns.Throughout the process the original L and r of a fibre (given in Table 2) has been kept constant.The data thus generated has been used for deriving the rest of our results.

Calculation of CD and CV
In system theory, a phase delay is the measure of how a certain frequency component of a signal is delayed by the system and therefore is more applicable for analysis of pure sinusoidal signals.On the other hand, the group delay gives the measure of the delay introduced by the system to the overall envelope of a signal composed of multiple frequency components.Since an AP is not a pure sinusoid but a composite signal, the time delay of its conduction through an INS must be estimated as a group Conduction Delay (CD) according to Eq. ( 14) given below: where θ(ω) is the frequency-dependent phase change introduced to an input signal by the system which could be directly calculated from the Bode Phase plot of the system TF.Therefore, for each TF created using the simulation flow mentioned above, we derived the CD directly from the corresponding Bode Phase plot by applying Eq. (10).Consequently, for an INS length L , CV could be calculated as CV = L/CD.

Figure 2 .
Figure 2. The variation of frequency-dependent gain characteristics of Aα11 fibre with changes of myelin turns ( M ).Filled circles denote f L in each case-note its shift toward the higher frequency range with increase in M and vice-versa.

Figure 3 .
Figure 3.The loci of f L with respect to variation in g and γ .All the PNS fibres show the same behaviour, and loci of their f L follow a linear curve.On the contrary, the CNS fibres show individual behaviour although still exhibit linear relation with g and γ given by Eq. (1).

Figure 4 .
Figure 4. Gain characteristics of the Aα11 fibre during compensation shows almost the same characteristics with f L coincident for all the curves.Solid circles show the f L for each TF calculated with r and L estimated using Eqs.(3a) and (4) for different values of M .Note that f L = 10 KHz is maintained even when M reaches to 2.

Figure 5 .
Figure5.The frequency-dependent CD (top) and corresponding CV (bottom) for Aα11 fibre in normal condition ( M = 400)-note that here we used angular frequency ω in the x-axis instead of linear frequency f as the computation is done according to Eq. (14); Section "Materials and methods, which requires representing the phase angle in the unit of radians.The filled circles denote ω L = 2πf L = 62832 rad/s corresponding to f L = 10 KHz.Note the value of CV below 1 kHz frequency ( ω ≤ 6283 rad/s) varies between 40 and 52 m/s-a value well within the experimentally found range since typically CV is measured within that frequency range.

Figure 6 .
Figure 6.Comparison of CD and CV during compensation.Although the CD (Top) at f L = 10 KHz( ω L = 62832rad/s ) remains constant, in both cases of M = 400 and 50 , the CV (Bottom) in the latter is much lower compared to the former.If measured within the range of 1 KHz ( ω ≤ 6283rad/s ), the CV will be even lower and consistent with experimentally measured range although this leads to a misinterpretation of the phenomenon.

Figure 7 .
Figure 7. Myelin geometry and corresponding electrical circuit model used for modelling.Top figure shows the cross sectional and longitudinal view of myelin structure with its geometrical parameters.In the bottom, R a , R m , R my and R p are the resistances of the axoplasm, axolemma, myelin sheath and periaxonal space respectively and C m and C my are the capacitances of the axolemma and myelin sheath.Inputs and outputs have been shown in complex s-plane.

Table 1 .
Electrical parameters used in this work.

Table 2 .
Geometrical data of the fibers under consideration.